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Abstract 

Inspired by a paper of R. W. Rosenthal, we investigate generalized Nash-equilibria of 
integer programming games. We show that generalized Nash-equilibria always exist and are 
related to an optimal solution of a so-called A^-fold integer program. This link allows us to 
establish some polynomial time complexity results about solving this optimization problem 
and its inverse counter-part. 



1 A class of IP games 

There are N players. For i = 1, . . . , TV, let A' e Z^'"^", b' e Z"^' , and u' e Z'^. Each player i selects 
a solution x* to the linear Diophantine system A'x — b',0 < x < u', x G Z", which is referred 
to as a pure strategy of player i. The costs to i for playing his strategy depend on the strategies 
chosen by the other players. A player is called satisfied if his selected strategy is cost- minimal 
given the strategies of all other players. A pure Nash- equilibrium is an assignment of strategies to 
players such that all players are satisfied simultaneously. Note that in the described situation each 
player can choose his strategy independently of the strategies of the other players. In contrast to 
this, one could impose a coupling constraint on the possible strategies for all players together, say 
B'x' < b° for matrices B' e Z™''", i = 1, . . . , iV, and for some b° e Z". Again, a player 
is called satisfied if his selected strategy is cost-minimal given the strategies of all other players. 
Clearly, his strategy has to obey the coupling constraint resulting from fixing the strategics of the 
other players. A generalized Nash- equilibrium is an assignment of strategies to players such that 
all players are satisfied simultaneously in this more general situation. Clearly, wc can model the 
unrestricted situtation by simply putting all to zero matrices. 

Let us be more specific now. Let there be given functions Cj : R+ that evaluate the cost of 

element j = 1, . . . , n. In our setting, we have costs Cj{y) if all the players together use y = J^iLi ^] 



1 



units of element j. In formulae, the costs assigned to player k are 

n / N \ n I N \ n / N \ n / N 

Y.-AY.y]] ~Y.-A E ^n=E^.U^+ E v)]-Y.-A E 

3=1 \i=i I j=i yi=i,i//c J j=i \ i=i,i=tk J j=i 

From now on, we assume that Cj is a monotonously increasing function for all j ~ 1, . . . ,n. The 
following result is an extension of a theorem by Rosenthal [5]. 

Theorem 1 All games of the above class with at least one choice (x^, . . . ,x^) of strategies that 
fulfill the coupling constraint < b*^ possess at least one generalized Nash- equilibrium. 

Proof. Consider the problem 

[E^i) ■■ A'x' = W,0 <x' <u\^B'^' <b°,x* E Z^l . (1) 
j=i \i=i ) 1=1 I 
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By assumption, this problem has at least one solution. Moreover, as for j = 1, . . . , n, the functions 
Cj are monotonously increasing, they are all bounded from below, and thus. Problem ([T|) has 
an optimal solution. We claim that any optimal solution to Problem ([T]) is a generalized Nash- 
equilibrium. 

Let (y^, . . . , y^) be an optimal solution to Problem For the purpose of deriving a contradiction, 
suppose that (y^, . . . ,y^) is not a generalized Nash-cquilibrium. Then there exists some player i, 
w.l.o.g. i = 1, and a strategy such that has lower costs 

n / ^ \ " / ^ \ 

3 = 1 \ 1=2 J 3 = 1 \i=2 J 

than the costs 

n / N \ n / N \ 



E^.N+E^} -E^. E; 



,y3 

j=l \ i=2 J 3 = 1 \i=2 I 

of strategy given fixed strategies y^, . . . , y^. We obtain 

n / ^ \ " ( ^ \ " / ^ \ " ( ^ \ 

E^.U^+^^0-^^4^^^r^'^4^^+^^0"^^4^^0' 

3 = 1 \ i=2 I 3 = 1 \i=2 I 3 = 1 \ 1=2 I 3 = 1 \i=2 J 

that is, 

Y.-A^^ + Y.y}]<T.-Ay^ + T.yU- 

3 = 1 \ 1=2 / j=l \ 1=2 / 

Thus, (z^, y^, . . . , y^) is a feasible solution to Problem ([1]) with strictly smaller costs than the 
solution (y^, y^, . . . , y^), contracting optimality. □ 

The following remark seems in order. Notice that the objective function in Problem(IT|) can be 
interpreted as the costs of a provider to account for realizing all the strategies. The costs assigned 
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to a player, however, are marginal costs. If all Cj, j = 1, . . . ,n, are nonlinear convex functions, 
the sum of costs for all players exceeds the costs for the provider, i.e., it generates a gain for the 
provider. The more players play the bigger this gain gets. This fact motivates the provider to have 
all his players being satisfied. Interestingly, each player is satisfied when they cooperate to minimize 
the costs of the provider. 



2 Complexity results for convex cost functions 

Let us assume in this section that all matrices are equal to A, that all matrices are equal 
to B, and that the Cj : M+ M+, j = are univariate convex functions. This latter 

condition is quite natural, since it simply means that the difference functions Cj{y + 1) — Cj{y) 
are monotonously increasing, that is, the cost for one extra unit of component j is monotonously 
increasing the more of component j is used. 

With the simplifications. Problem ([T]) has a nice combinatorial structure: 
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This structure allows us to prove a complexity result about the computation of the generalized 
Nash-equilibrium. However, in order to prove this, we need to introduce the notion of the Graver 
basis of a matrix together with a few important facts about these bases and about A^-fold integer 
programs. 

For D G Z'^^" let G{D) denote the Graver basis of D. This Graver basis consists exactly of those 
nonzero elements z in ker(A) fl Z" = {x : Ax = 0} fl Z" that cannot be written as a non-trivial 
sum z = u + V where u, v S ker{A) n Z" and where u and v have the same sign-pattern, that is, 
UjVj > for all j = 1, . . . ,7i. For any matrix D, the Graver basis 0(0) is a finite set [5], and it 
has the nice property that it provides an optimality certificate for the family of separable convex 
integer minimization problems min |X]j=i fji^j) '■ = d,0<x<u, x£ | for any choice of 

convex functions /j : M ^ M, j = 1, . . . , n, and for any choice of d G Z'* [11 [5]. That is, for any 
choice of the input data and for any feasible solution x" of the problem, x*' is optimal if and only 
if there is no g G GiD) such that x° + g is feasible and J27=j fii^^ + 9j) < 51]"=^ fji^j)- 

Having a finite optimality certificate available, one can augment any given feasible solution of the 
problem to an optimal solution by iteratively walking along directions provided by the optimality 
certificate. If a simple greedy augmentation is followed, only a polynomial number (in the encoding 
length of the input data) of Graver basis augmentation steps is needed to construct an initial feasible 
solution and to augment it to optimality [U Thus, if we had a Graver basis G{D) that is of 
polynomial size in (D) , this augmentation approach provides a polynomial time algorithm to solve 

min 1^;^, /j (xj ) : L'x = d, < X < u, X G Z!^ } . 
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Using results of [H [7] , it was shown in [T] that for fixed matrices A and B (of appropriate dimen- 
sions) the sizes and the encoding lengths of the Graver bases of the family of matrices 



[A,B] 
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increase only polynomially in A^. Thus, if we choose D = B]^^^ the corresponding separable 
convex minimization problem, a so-called A'^-fold integer program, can be solved in polynomial 
time for fixed matrices A and i? [U [3] . The following lemma shows that for fixed matrices A and 
B also the Graver bases of Problem ^ increase only polynomially in N . And thus Problem ^ 
can be solved in polynomial time as claimed. 



Lemma 2 Let A G Z''^" and B G j/^Lxn fixed matrices, and let 
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Then the sizes of the Graver bases of [A, B] , N = 1,2, . . and their encoding lengths increase 
only polynomially in N . 



Proof. Consider the matrix 
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As this is the matrix of an iV-fold IP composed of block matrices ( A o o ) and o o /„ , the sizes 

\ B / 

of the Graver bases of C'^\ iV = 1,2,..., and their encoding lengths increase only polynomially 
in A^, given that A, B, and n are constant. 

{N) 



is composed out of columns of C^^^ and thus any Graver basis element of 
filled with zeros for the additional columns of (7^^-* defines a Graver basis element of 



The matrix [A, B] 

C(^). Therefore, also the sizes of the Graver bases of [A, i?] iV = 1, 2, . . ., and their encoding 
lengths increase only polynomially in TV, given that A, B, and n are constant. □ 
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Theorem 3 Assume that A G l^dxn g Z™'^" are fixed and assume that all cj : M+ R+, 

j ~ l,...,n, are univariate convex functions given by evaluation oracles. Moreover, let M be a 
bound for Cj (^J2iLi f'^^ 3 '^'^^ /'^'^ feasible selections of strategies for all players. Then 
Problem can be solved, and hence a generalized Nash- equilibrium can be computed, in time that 
is polynomial in N , in {M, , b"'^, . . . , b^, U"'^, . . . , U"'^), and in the number of calls of evaluation 
oracles. 



Proof. From Lemma HJ we know that the Graver basis of the problem matrix [A, B] has a 
polynomial-size encoding length in the encoding length of the input data. The result now follows 
from Theorem 7 in [3], which bounds the number of augmentations so as to reach an optimal 
solution of Problem □ 



3 Inverse optimization problem 

A (generalized) Nash-equilibrium satisfies all players simultaneously. It requires that cost functions 
Cj, j = I, . . . ,n, are given. Conversely, one may ask for finding cost functions cj, j ~ 1, . . . , n, that 
turn a given set {x^,...,x^} of strategies into a Nash-equilibrium. Such inverse optimization 
problems are fundamental for a variety of applications, such as quality control in production 
systems, portfolio optimization, or production capacity planning, see [S] for details and references. 

The problem in this generality is hopeless unless some structure on the cost functions Cj, j = 
1, . . . ,n, can be assumed. We study the situation that each Cj has the form Cj{y) ~ Xjfj{y) for 
some fixed given convex function fj : IR+ IR+ and some e IR+ that can be adjusted. 

We can abstract from this setting slightly and define a general inverse integer optimization problem 
as follows. 



Lemma 4 Consider the problem (HOP): Given a polyhedron P = {x : Dx = d, < x < u} C M'\ 
X* G P n Z" and convex functions fj : R — > M. Do there exist Xj £ j = 1, . . . ,n, not all 
simultaneously zero, such that x* minimizes /(x) :~ X]j=i ^jfj^-'^j) over all x P f] Z". 

(a) (HOP) IS m coNP. 

(b) (HOP) is solvable in time polynomial in {Q{D)), where Q{D) denotes the Graver basis of D. 

Proof. By Theorem 3 in [3], x* minimizes /(x) = X]j=i ^jfji^j) o'^*^'" all x € P n Z" if and if 
only if /(x* -f g) > ./(x*) for all g g GiD) such that x* + g e P. Clearly, this can be checked in 
polynomial time in {Q{D)), showing (b). 

To show (a), observe that this criterion is equivalent to asking for a solution to the following system 
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in variables Ai, . . . , A^: 

n n 

^ \jfj{x* + gj) > Xjfj{x*) for all g £ g{D) such that x* + g e P, 
j=i i=i 
Ai, . . . , A„ > 0, 

n 

^A, > 0, 

i=i 

or equivalcntly, letting if = {g S 5(1?) such that x* + g G P}, 

n 

^ foraUgeif, 
Ai, . . . , A„ > 0, 

n 

^A, > 0. 

If this system has no solution, then there exists a subsystem 

n 

E[/^(^J*+3^)-/^(^p]^^- ^ foraUgei/', 
Ai, . . . , A„ > 0, 

n 

^A, > 0, 

for H' C H oi at most < n + 1 inequalities having no solution. (This is a consequence of 
Caratheodory's theorem and the Farkas-Lemma, see p. 94 in [5].) 

Now we can apply standard linear programming technology: From the Farkas-Lemma, it follows 
that the latter system is not solvable if and only if there exists a vector v = {vg)g(zH' & such 
that 

E (^J + ff^- ) - (^J )] < for all j = 1 , . . . , n. 

seH' 

This is a polynomially sized certificate for the "No. "-answer to (HOP), showing (a). □ 

Notice that Lemma|4](b) implies that if a Graver basis happens to be of polynomial size in the en- 
coding length of the input data {D, u) of the integer program, then the corresponding optimization 
problem (HOP) can be solved in polynomial time in (Z?, u). 

Returning to Problem we want to solve its inverse counter-part. Note that again we assume 
that each Cj has the form Cj{y) = ^jfjiu) for some fixed given convex function fj : M_|- —> M_|- and 
some Xj G K+ that can be adjusted. 

Corollary 5 Let the matrices A and B be fixed and assume that all fj : R+ —> R+, j ~ 1, . . . ,n, 
are univariate convex functions given by evaluation oracles. Moreover, let M be a bound for 
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^jj '^^^ '^'^'^ /'^^ feasible selections of strategies for all players. Then the in- 
verse integer optimization problem corresponding to Problem 0j is solvable in time polynomial in 
N , in (M, X*, b", b"'^, . . . , b"'^, U"'^, . . . , u^), and in the number of calls of evaluation oracles. 

Proof. From Lcmnia[21 we know that the Graver basis of the problem matrix [A, B] of Problem 
([2]) is of polynomial size in A^. The result thus follows from Theorem [Jb) . □ 



4 Extension of results 

Let us call the pair of matrices (j4\ i?') the type of player i. While in Problem ([T]) all players could 
have a different type, we assumed in Problem ([2]) that all players were of the same type {A,B). 
Let us assume now that there is a constant number of possible player types represented 
by pairs of matrices {A^,B^), i = 1,. ..,t, as in Problem ([1]). Despite the more general problem, 
the complexity results of the last two sections still hold. 

Theorem 6 Assume that there is a constant number t of possible player types represented by pairs 
(A\ B''), i — 1, . . . , of fixed matrices A^ and B^ , i ~ For i = 1, . . . , N , select one of 

these pairs B"^), i = 1, . . . , i, as the type (A*, B^) of player i. Finally, let M be a bound for 

Cj (^J2iLi fof 0,11 j o,nd for all feasible selections of strategies for all players. 



(a) Problem {Ip can be solved in time polynomial in N, in (M, b'^,b^, 
in the number of calls of evaluation oracles. 



), and 



(b) Assume that each Cj has the form Cj{y) = ^jfj{y) for some fixed convex function fj : M-|- —^ 
R+ given by evaluation oracles and some Xj E M+ that can be adjusted. Then the inverse 
integer optimization problem corresponding to Problem (QP is solvable in time polynomial in 
N, m (M,x*,b°,b\...,b^,ui,...,u^ ), and in the number of calls of evaluation oracles. 



Proof. As in the proofs above, it suffices to show that under our assumptions the sizes of the 
Graver bases of the problem matrices of Problem ^ increase only polynomially with N. This can 
be seen by considering the matrices 
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which are matrices of iV-fold integer programs. Thus, as these A''-fold matrices are built from fixed 
blocks, the sizes and encoding lengths their Graver bases increase only polynomially with N. As 
the problem matrix of Problem ([1]) can be obtained from the matrix by deleting columns and 
then zero rows, the sizes and the encoding lengths of the Graver bases of the problem matrices of 
Problem ([T|) increase only polynomially with N, and thus, (a) and (b) both follow. □ 
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